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The gravitational waveforms of a chaotic system will exhibit sensitive dependence on initial con- 
ditions. The waveforms of nearby orbits decohere on a timescale fixed by the largest Lyapunov 
exponent of the orbit. The loss of coherence has important observational consequences for systems 
where the Lyapunov timescale is short compared to the chirp timescale. Detectors that rely on 
matched filtering techniques will be unable to detect gravitational waves from these systems. 
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The evolution of two compact objects under their mu- 
tual gravitational attraction is a major unsolved problem 
in general relativity. A vast array of analytic and numeri- 
cal approximations have been used to attack the problem, 
but much remains to be understood. Amongst the many 
outstanding questions is the degree to which the dynam- 
ics may be chaotic §, § |, § or effectively chaotic § |. 
For example, when one or more of the masses is spin- 
ning it has been shown in the test-particle [^| and post- 
Newtonian [0 limits that certain classes of orbits are 
chaotic, at least when radiation reaction is turned offQ. 

Here we consider the physical and observational con- 
sequences wrought by chaotic behaviour in the orbital 
dynamics. A key feature of chaotic systems is their sen- 
sitive dependence on initial conditions - the so-called 
"butterfly effect". It is intuitively obvious that sensi- 
tive dependence in the dynamics will be reflected in the 
gravitational waveforms, but the precise connection has 
not been established until now. We find that the wave- 
forms of nearby orbits decohere on a timescale fixed by 
the largest Lyapunov exponent. The observational conse- 
quences of this result depend on two timescales: the chirp 
timescale T c ; and the Lyapunov timescale T\. A third 
timescale of interest is the gravitational wave period T w . 
The wave period and chirp timescale for a system with 
reduced mass /i, total mass M and orbital separation R 
are given by 



T„, ~ ttM 



3/2 



and T r 



256 n 



In the limit that one of the masses is very much smaller 
than its companion we have (i <C M and T c 3> T w . The 
Lyapunov timescale has to be calculated on a case-by- 
case basis, but values as short as T\ ~ T w are possi- 
ble for very unstable orbits. Orbits whose Lyapunov 
timescale is short compared to the chirp timescale will 
produces highly unpredictable waveforms. The number 
of templates required to detect these waveforms is expo- 
nentially large, making them impossible to detect ||, 0] 
with the Laser Interferometer Gravitational Observatory 
(LIGO). 

The connection between orbital instability and wave- 
form decoherence is easily established. It relies on the 
fact that the gravitational waveform can be ex- 

pressed in terms of the phase space coordinates of the 
system, so the divergence of nearby trajectories in phase 
space is reflected in the divergence of nearby waveforms. 



The phase space coordinates we have in mind are the 
position and momentum of the masses. These coordi- 
nates are well defined in the test particle, post-Newtonian 
and post- Minkowski approximations, but their meaning 
is less clear in the full theory. Indeed, the full evolution 
equations are a set of coupled, non-linear partial differ- 
ential equations, so the real issue is gravitational turbu- 
lence, not chaos. We will avoid this complication by as- 
suming that the evolution equations can be approximated 
by some set of non-linear ordinary differential equations 
for the phase space variables Xi(t). Here t is taken to be 
the time measured by a distant observer, but the choice 
is unimportant. The evolution equations can be written: 



dt 



(1) 



Linearizing about a solution to this equation yields an 
equation for the perturbation: 



d6Xi{t) 
dt 



= K ij (t)SX j (t). 



where 



Kij{t) = 



dHi 



dXj 



(2) 



(3) 



is the infinitesimal evolution matrix. The solution to the 
linearized equations of motion can be expressed in terms 
of the evolution matrix Ly(t): 



dX i (t) = L ij (t)dX j (0). 



(4) 



The eigenvalues and eigenvectors of K+j (t) contain infor- 
mation about the local stability of the phase space trajec- 
tory. In a fixed coordinate system the eigenvectors may 
vary with time. To avoid this, a non-constant basis can 
be defined that stays aligned with the eigenvectors^). 
This is done by parallel transporting the basis vectors 
along the trajectory Xi(t). Barring degeneracies, the in- 
finitesimal evolution matrix is diagonal in the eigenvector 
basis: Kij(t) — Sij\i(t) (no summation on the i), where 
the Ai(t) are the eigenvalues of Kijit). The evolution 
matrix is also diagonal in the eigenvalue basis and has 
components 



Lij(t) = S tj exp 



\i(t')dt' 



(5) 
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The principal Lyapunov exponent for the trajectory 
Xi(t) is defined: 



1 



(6) 



The multiplicative ergodic theorems establish the exis- 
tence of this limit for a large range of situations. Tra- 
jectories with positive Lyapunov exponents are unstable 
against small perturbations. When the unstable orbits 
have non-zero measure in phase space, the system is said 
to be chaotic. A positive principal Lyapunov exponent 
will dominate the long term dynamics. Setting Ai to be 
the principle eigenvalue, the asymptotic form of the evo- 
lution matrix can be written: 



Lij(t) 



/(*) o o 






(7) 



where \f(t)\ = 1. If X x {t) is real, then f(t) = 1. 

A gravitational wave can be decomposed into two po- 
larizations with amplitudes h (t) and h x (t). A gravi- 
tational wave detector responds to the wave according 
to 



»(t) = h+(t)F + +h x (t)F x , 



(8) 



where F + and F x are the antenna patterns for each po- 
larization. The antenna patterns for LIGO and LISA can 
be found in Refs.||] and |l(| respectively. The gravita- 
tional wave amplitude h(t) measured at the detector can 
be expressed in terms of the phase space variables Xi (t) . 
Nearby orbits have gravitational waveforms that differ by 



Sh(t) 



dh(X m ) 
dX t 

dh(X m ) 
dX, 



(t)SX t (t) 

WLvMSXjiO). (9) 



Employing an eigenvector basis and using the asymptotic 
form for the evolution matrix yields 



where 



Sh{t) 



dh(X v 



dX x 



-(t)f(t) 5X^0) 



(10) 



(11) 



is an oscillatory factor with bounded amplitude. The 
divergence of the waveforms can be better expressed in 
terms of their relative phase. Writing the amplitude of 
the reference trajectory as and defining $o(t) to be 

the phase for which h(&o) = 0, we find the relative phase 
is given by 



5®(t) 



&'(#„(«)) 



(12) 



The quantity 



A/j = lim lim In 



6$(t) 



5<Z>(0) 



(13) 



is equal to the principal Lyapunov exponent of the tra- 
jectory. In other words, the waveforms of nearby trajec- 
tories will decohere exponentially in time, on a timescale 
set by T\ = A" 1 . 

A couple of comments are in order. Firstly, the limit 
t — > oo used in equations (||) and (13) will return 
A = Afj = if we include the effects of radiation re- 
action. This reflects the fact that the dynamics is not 
chaotic in a formal sense Q. The limit t — > oo has to 
be replaced by the limit t -> T, where T\ < T < T c , 
This is done by replacing the Lyapunov exponent by a 
finite-time or local Lyapunov exponent and making a 
similar change to (|l|). Clearly, the whole notion of wave- 
form decoherence only makes sense if the Lyapunov time 
scale T\ is very much shorter than the chirp timescale 
T c . It remains to be seen if this condition is satisfied 
for any realistic systems. The second comment concerns 
the gauge non-invariance of the Lyapunov exponents and 
the choice of time variable. Lyapunov exponents gained a 
bad reputation when they were used to study the Bianchi 
IX dynamics. Different time slices yielded principal Lya- 
punov exponents that were positive or zero, making it 
difficult to decide if the dynamics was chaotic. However, 
this difficulty can be avoided by making relative, rather 
than absolute comparisons. For example, the ratio T\/T c 
is the same in any coordinate system, so the choice of 
time variable is irrelevant. The confusion surrounding 
the Bianchi IX dynamics could have been avoided if the 
Lyapunov timescale had been compared to the average 
time between bounces [[l2| . 

We can illustrate the connection between orbital in- 
stability and waveform decoherence with a simple exam- 
ple. Consider the Lagrangian for a test particle in the 
Schwarzschild spacetime (in units where G = c = M = 1) 




(14) 



The system is completely integrable (non-chaotic) as it 
has four generalized coordinates t, r, 9, <p and four con- 
stants of motion, the mass [A, energy E, z-component 
of angular momentum L z and total angular momentum 
squared L 2 . Despite the lack of chaos, the dynamics does 
admit isolated unstable orbits that serve to illustrate the 
connection between orbital instability and waveform de- 
coherence. Restricting our attention to equatorial orbits 
with fixed energy E allows us to make the replacement 



2 d 



dX 



rE dt 



(15) 



in the equations of motion: 
E 3 r 

Pr -- 



Ep 2 r 



(r-2) 3 r(r-2) r 2 (r-2) 
P<t> = 
r = i?j) r 
• = E P4> 
9 r(r - 2) ' 

Our reference trajectories are circular orbits with 

Pr = 

= L z = L 

r = r± = ^L(L ± V ' L 2 - 12) 
= ujt. 



(16) 



(17) 



Perturbing the equations of motion about the reference 
trajectory yields the infinitesimal evolution matrix 



Kij = 





E 




2EL 

r 2 (r-2) 




E 



r(r-2) 

The eigenvalues are: 



2£ J (r+l 

(r-2) 4 r 3 (r-2) a 




2EL(r-l) 
r^r-2) 2 



(18) 



Ai — — A2 — 
A 3 = A 4 = 



y/2E 2 r 3 {r + 1) - L 2 (3r - 4)(r - 2) 2 



r 2 (r- 2) 2 



(19) 



The well known innermost stable circular orbit at r = 6 
has L = 2V3, E = ^/8/9 and Ai = A 2 = 0. Circular 
orbits with larger values of L come in two varieties, those 
with r = r + and those with r = r_ , e.g. setting E — 1 
and L = 4 yields orbits with r + = 12 and r_ = 4. The 
outer orbit is stable, A x = — A 2 = 1^6/108, while the 
inner orbit is unstable, Ai = — A 2 = 1/2^/2. The unstable 
orbit at r — 4 has eigenvectors 



ei 



e 2 = 



1 2 1 

372'°'3' _ V2 

3V2 ' 3' V2 



e 3 = e 4 = (1,0,0,0). 



(20) 



Since the zero eigenvalues have degenerate eigenvectors, 

can not be diagonalized, but it is enough to put 
into Jordan normal form with the transformation matrix 



Pij = 



TVS 


7>/2 


24 


24 








7 


7 


6 


6 




1 

-2 



(21) 



In terms of the transformed phase space coordinates the 
infinitesimal evolution matrix has components 



K, 



and the evolution matrix has components 



1 

2V2 














1_ 

















1 















(22) 



Lij — 



a t/2y/2 







,-t/2V2 



1 t 











(23) 



Taking the limit defined in (|g) yields the principal Lya- 
punov exponent A = l/2v2. The Lyapunov timescale 
for this orbit, T\ = A -1 = 2^/2, is actually shorter than 
the gravitational wave period T w = 27r. 

Next we need an expression for the gravitational wave 
amplitude in terms of the phase space variables. The cor- 
rect approach would be to solve the Teukolsky equation 
using the method described by Hughes Q, but no ana- 
lytic solutions exist in the strong field limit. Instead we 
will use the quadrupole approximation, knowing full well 
that it is being applied outside of its domain of validity. 
Since our goal is to illustrate a general phenomenon - not 
to produce accurate templates for LIGO data analysis - a 
qualitative description of the waveform's dependence on 
the phase space variables will suffice. 

The gravitational wave amplitude can be written as 
h^uit) — h + {t)e+ v + h x (t)e* v , where the polarization 
tensors have non-zero components 

e+ = 1 

1 



"yy 



■xy 



-yx 



1 . 



(24) 



A detector situated at r = R on the z axis will encounter 
a wave with plus polarization given by 



h+(t) = 



2fiE 2 
R 



^P<j>Pr{r - 3) 



sin 2(f> - 



E z r 



2J2 



P%r + 2) p 2 (r-3) 
r(r - 2) 2 r-2 



cos 2(f) 



(25) 



and similarly for the cross polarization. Specializing to 
the orbit with E = 1, L = 4 and r = 4 we have 



h+ 
dh+ 



8/i 



R 
2u . 



cos 20 



4 



so that 



dh+ 


3/i 


dp<t> 




dh+ 


17^ 


dr 




2R 


dh+ 




deb 


1T ! 



\(b. (26) 
Putting everything together in equation (||) yields 

5h+(t) ~ ^e* /2 ^(119cos2</>- 175\/2sin2</>) 

x(3<5r(0) +6(5p (O) + 6V2<5p r (0)), (27) 



<y$(t) 



,t/2A/2, 



<5$(0) . 



(28) 



As promised, the phase difference grows exponentially on 
a timescale equal to T\ = 2^/2. 
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